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Dividing by mixed-events for acceptance correction is wrong 
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We point out that the mixed-event method for two-particle acceptance correction, widely used in 
particle correlation measurements at RHIC and LHC, is wrong. The correct acceptance should be 
the convolution of two single-particle efficiency x acceptance functions. The error of the mixed-event 
method, which guarantees a uniform A77 two-particle combinatorial density, is, however, small in 
correlation analyses where the two particles are integrated over an extended pseudo-rapidity rj range. 
With one particle fixed in 77 and the right acceptance correction, the background-subtracted corre- 
lated pair density may reveal not only a short-range but also a long-range Arj dependence. This has 
important physics implication, and may provide crucial information to disentangle physics mecha- 
nisms for the recently observed long-range ridge correlation in asymmetric proton-lead collisions at 
the LHC. 

PACS numbers: 25.75.-q, 25.75.Dw 
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I. TWO-PARTICLE ACCEPTANCE 

Two-particle correlations are a valuable tool to study 
heavy-ion collisions [TH5]. Correlation functions are of- 
ten formed by particle pair density in real events divided 
by that from mixed-events, where the two particles are 
taken from different events. With proper normalization 
the deviation of the correlation function from unity re- 
veals correlations between two particles. This technique 
applies to situations where all particle pairs are corre- 
lated, for example, in Hanbury-Brown and Twiss (HBT) 
interferometry [3] and anisotropic flow correlations [S]. 

One unique case of correlations is jet angular correla- 
tions. The study of jet correlations has provided wealth 
of information about relativistic heavy-ion collisions |T]- 
|3]. The object of interest is a cluster of particles-not 
all particles in the event -that are correlated due to their 
common origin of parton fragmentation. The interest is 
often to find number of correlated particles and their an- 
gular distributions. The correlations are often formed 
in terms of the particle pair azimuthal angle difference, 
A(J3 and pseudo-rapidity difference, Arj. The mixed-event 
technique is widely used, not as the uncorrclatcd baseline 
as in HBT or anisotropic flow analysis, but to correct for 
two-particle acceptance. For example, in mid-rapidity 
region, the mixed-event two-particle density is approxi- 
mately triangular in A77, and is used for Arj acceptance 
correction after normalized to 100% at Ar] = 0. This 
mixed-event method for two-particle acceptance correc- 
tion is, however, wrong. The correct two-particle ac- 
ceptance should be the convolution of two single-particle 
acceptances (or efficiency functions). It has nothing to 
do with the measured single-particle densities, but only 
depends on detector efficiencies. The mixed-event two- 
particle density, on the other hand, depends on both the 
detector efficiencies and the true single-particle densities. 

In the widely used mixed-event acceptance correc- 
tion, the single-particle distributions are mistaken as de- 
tection efficiencies, implying that the true distributions 
(e.g. the pseudo-rapidity density dN/drf) are always uni- 



form. By using mixed-event technique, the "corrected" 
two-particle correlation signal in Arj is guaranteed to 
be uniform except regions of correlation signals. The 
truth may, of course, not be uniform as the single-particle 
dN/dr] at mid-rapidity is not strictly uniform. This, how- 
ever, does not make a big error for mid-rapidity particles 
in symmetric collisions, where particle density is nearly 
uniform in ry 0. 

The mixed-event acceptance correction could be fur- 
ther from truth for asymmetric collisions, such as proton- 
lead (p-Pb) collisions, where the single-particle dN/dr] is 
nonuniform even at mid-rapidity ^. The corrected two- 
particle correlations will, again, be uniform by construc- 
tion except the A77 regions of any correlation signal. But 
this clearly is wrong. It is easy to see in the following sim- 
ple example. Suppose in a jet-correlation study in A77, 
the high-pT trigger particle is fixed at 77 = and all other 
particles are paired with the trigger to form angular cor- 
relation function. There will be a peak at Arj = due to 
jet correlations and the underlying background will have 
the same shape as the underlying event single-particle 
dN/dr] distribution. The two-particle acceptance correc- 
tion to be applied should be the single-particle efficiency 
as function of 77 (because the trigger particle is always at 
T] — 0). After correction, the signal should be the real 
jet signal on top of the real background which, in this 
case, is the true single-particle density which may not be 
uniform in r/. If the mixed-event A?7 distribution (which 
is the single-particle rj distribution in this simple case) 
is used as the acceptance correction, then the corrected 
pair density will be completely flat in A77 except for re- 
gions of correlation signal. Clearly, this does not reflect 
the true physics condition and will be wrong. 

In real data analysis, often all pairs of trigger and 
associated particles are used. In asymmetric collisions, 
neither the trigger nor the associated particle rj density 
may be uniform. Averaging over all trigger and associ- 
ated particles, the non-uniformities in the single-particle 
dN/dr] distributions, however, become a second-order ef- 
fect in the two-particle density distribution of dN/dArj. 



To see this, we take the simple example of a linear, mea- 
sured single-particle density. 
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where ±77™ are the acceptance limits. The combinatorial 
two-particle density distribution will be 



dN 
dArj 



m 






(27;„-|A77|) 



l + k 

niax( — ry^i , — 7/^11 — Ar/) \ '/' 

1 



1 + k — 1 5{i-j2 - til - A7])dr]idr]2 



1-f -fcM 2-2 
6 \ -q, 



r 




(2) 



The effect of non-uniformity is quadratic in k. When 
k « 1 which is typically the case for mid-rapidity re- 
gion, the two-particle density is approximately triangular 
in Ar/. Note that the two-particle density integrated over 
trigger and associated particle 77 in Eq. ^ is symmetric 
with respect to A77 = even though those for fixed trig- 
ger r] are all asymmetric. This is because, taking fc > as 
example, the larger associated particle density at positive 
A77 is weighted by the relatively fewer trigger particles at 
negative 77, and vise versa. 

In most A77-A0 correlation analyses, the mixed-event 
technique is used for two-particle acceptance correction 
and, as we discussed, the final correlated yield is uni- 
form by construction. If the correct acceptance is used 
from convolution of single-particle efficiencies, the result- 
ing two-particle density does not deviate much from a 
uniform distribution, as we have shown above. For ex- 
ample, for a perfect single-particle efficiency, the real two- 
particle acceptance is a triangle. The relative error the 
mixed-event correction makes by using (the properly nor- 
malized) Eq. (pi) is 
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For relatively small k, the error due to mixed-event two- 
particle acceptance correction is small. For fc « 1/4 as 
measured in central 20% d-|-Au collisions at RHIC by 
PHOBOS over a hypothetical acceptance range of 
— 2 < 77 < 2 (i.e. rjm = 2), Eq. ([3| gives ^ 3% error 
at A?] — rjrn and ^ 8% error at Arj = 277^- This is 
illustrated in Fig. [T] 



II. PHYSICS IMPLICATION 

We have so far focused on a technical point, that two- 
particle acceptance correction in A77 should not be ob- 
tained from mixed-events, but rather from convolution of 
two single-particle efficiency x acceptance functions. Now 
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FIG. 1: Relative error given by Eq. ([3| that the mixed- 
event acceptance (Cmix) makes relative to the true two- 
particle acceptance (Ctriangie) for a perfect detector (100% 
efficiency X acceptance uniform in 77). 



we want to turn to an important physics implication of 
this point. 

Surprisingly, a strong large-A77 (A77 > 2) near-side 
(Acf) « 0) correlation is observed in p-Pb collisions above 
a uniform ZYAM (zero yield at minimum ^) background 
at the LHC [5HII]- This is called the "ridge" following 
the observation in heavy-ion collisions where an elliptic 
flow modulated background is subtracted [T2HT4J. The 
strength of the ridge in p-Pb is as strong as that in heavy- 
ion collisions, where the ridge is considered to be primar- 
ily a consequence of anisotropic triangular flow. As we 
discussed, the observed ridge in p-Pb collisions is uniform 
in A77 by construction of the mixed-event technique, and 
will likely remain uniform even using the correct two- 
particle acceptance correction because the trigger and 
associated particles in a wide 77 range are averaged. How- 
ever, the extra feature of the nonuniform single-particle 



1] distribution in p-Pb collisions may be essential to un- 
ravel the underlying physics mechanisms for the ridge 
correlations. In order for any possible A?7-dependence of 
the ridge to be observable, one needs to fix the trigger 
particle rjtrig within a narrow bin and study correlations 
of associated particles, dN/dAr] « dN/dr]assoc- 

There are two leading theoretical models for the 
physics mechanism of the ridge in small systems. One 
is hydrodynamics where the initial geometry anisotropy 
is converted into final-state anisotropic (dominantly el- 
liptic) particle distributions by hydrodynamical evolu- 
tion [15 . The anisotropic particle distributions result 
in an enhanced two-particle density at Acj) — 0, mani- 
fested as the ridge. In this picture, the ridge strength is 
proportional to the underlying background pair density: 
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Because the underlying event baseline dN (Arj) / dArj de- 
pends on A77, a measurement of the signal d^ N / dAr]dA4) 
will be informative about the nature of the ridge. Of 
course, the anisotropic harmonic Vn (product of the two 
single-particle azimuthal anisotropics, w„) may also de- 
pend on A77. Thus, a measurement of d'^N/dArjdAi] may 
not uniquely confirms or refutes the hydrodynamic expla- 
nation, but should provide an important extra informa- 
tion. In this respect, it would be crucial to have predic- 
tions of the 77 dependences of V2 and W3 from hydrody- 
namical calculations. 

The second model is the color glass condensate 
(CGC) [T5]. The CGC framework gives particular predic- 
tions of the two-gluon production process as a function 
of Arj. A calculation for Au-|-Au collisions is given in 
Ref. 17 which attributes the ridge to a net effect of the 
CGC enhanced two-gluon production and the strong col- 
lective radial flow. Calculations of the Arj dependence 
of the ridge in p-Pb collisions are not yet available but 
should be extremely valuable. 

Presumably, different physics mechanisms would yield 
different At^ dependences of the ridge correlations, as well 
as different energy dependences. Measurements of the 
ridge as a function of A77 at both the LHC and RHIC 
should, therefore, put stringent constraints on models. 

Experimentally, the two-particle acceptance should be 
taken from the convolution of two single-particle efficien- 
cies. For trigger particles in a fixed narrow 77 range, 
the two-particle acceptance is approximately the single- 



particle efficiency X acceptance for associated particles. 
One obtains the near-side raw yield vs. A77 and sub- 
tracts the combinatorial background by, for example, the 
ZYAM procedure [S]. The ZYAM magnitude as a func- 
tion of Ai] may be treated as the underlying background 
particle pair density. One then examines whether the 
background-subtracted signal is uniform in Ar/, propor- 
tional to the background pair A77 density in the ZYAM 
region, or of any other shape. 

To examine the ridge yield relative to the combina- 
torial background, one may not even need to do two- 
particle acceptance correction. One can simply form the 
ratio of the pair density from real events to that from 
mixed-events (i.e. the original correlation analysis), but 
with fixed narrow rj bins for trigger particles, and ex- 
amine the shape of the correlation function in the ridge 
region of large Arj. One may not even need mixed-events, 
but treating the pair density in the ZYAM A0 region as 
the combinatorial background, and examine the ratio of 
the pair density in the ridge region at small A(j) to that 
in the ZYAM region, both of a function of Arj. 

The idea of Ary dependence of the ridge can be ap- 
plied to asymmetric heavy-ion collisions, for example, the 
available Cu-f Au collision data from RHIC. 

In summary, we point out the misconception in two- 
particle acceptance correction by the mixed-event tech- 
nique. The mixed-event acceptance correction is, in 
principle, wrong and guarantees a uniform A77 pair 
density (correlation) signal. The correct two-particle 
acceptance should be the convolution of two single- 
particle efficiency X acceptance functions. The effect of 
the improper mixed-event correction is, however, small in 
present experimental data where both the trigger and as- 
sociated particles are integrated over extended ij ranges. 
With a fixed narrow 1] range for trigger particles and the 
proper acceptance correction, the correlation signal may 
reveal a A77 dependence, especially for asymmetric colli- 
sion systems such as p-Pb collisions recently available at 
the LHC. A A77 dependence, if present, could be crucial 
to discriminate physics models for the observed ridge. 
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